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Due to increasing research being conducted at ICASE in the field of fluid mechanics, 
future ICASE reports in this area of research will be printed with a green cover. Applied 
and numerical mathematics reports will have the familiar blue cover, while computer science 
reports will have yellow covers. In all other aspects the reports will remain the same; in 
particular, they will continue to be submitted to the appropriate journals or conferences for 
formal publication. 
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and 
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ABSTRACT 

Simplified equations for slender jets are derived for a circular jet of one fluid flowing into 
an ambient second fluid, the flow being confined in a circular tank. Inviscid flows are studied 
which include both surface tension effects and Kelvin-Helmholtz instability. For slender jets 
a coupled nonlinear system of equations is found for the jet shape and the axial velocity 
jump across it. The equations can break down after a finite time and similarity solutions 
are constructed, and studied analytically and numerically. The break-ups found pertain to 
the jet pinching after a finite time, without violation of the slender jet ansatz. The system 
is conservative and admissible singular solutions are those which conserve the total energy, 
mass and momentum. Such solutions are constructed analytically and numerically, and in 
the case of vortex sheets with no surface tension certain solutions are given in closed form. 
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1 Introduction 

Rayleigh [1] first studied the break-up of liquid jets in air. The underlying instability mech- 
anism is a consequence of surface tension, and Rayleigh predicted that all waves longer than 
the jet diameter are linearly unstable and used the result to conjecture the break-down of 
the jet into drops. The physics of the instability can be understood in purely energetic 
terms since for a given stationary volume of fluid the surface energy is a minimum when 
the shape is spherical. A liquid thread, therefore, will tend to break up into drops. An 
equivalent and more quantitative feature of the instability is found also by consideration of 
the Young-Laplace equation which relates the pressure jump across the jet interface to its 
curvature. In cylindrical geometries we have two principal radii of curvature, one in a plane, 
containing the axis of the jet and one in a plane perpendicular to the axis. If we suppose 
that a perturbation causes the jet radius to decrease slightly, then the radius of curvature 
in the cross-sectional plane decreases and the pressure in the vicinity of the perturbation 
increases. This means that fluid will tend to leave this high pressure region to flow to lower 
pressure positions, which by mass conservation decreases the jet radius further thus enhanc- 
ing the interfacial instability. This is a physical description of the mechanism of pinching 
or break-up. Clearly the phenomenon is nonlinear and even though linear theory does well 
in predicting dominant wavelengths in natural transition, nonlinear analyses are required to 
follow the jet to break-up for a range of dominant length scales whether self-sustained or 
externally imposed through initial conditions for example. 

The physical picture sketched above is fundamental and appears in both inviscid as well 
as viscous flows. In this work we concentrate on inviscid irrotational flows of liquid jets. A 
review article has been written by Bogy [2]. Of relevance is also the work of Chandrasekhar 

[3] and some experimental pictures of jet break-up can be found in Donnelly and Glaberson 

[4] (see also Drazin and Reid [5]). A theoretical study of the break-up of an inviscid incom- 
pressible liquid jet in air with surface tension has been carried out by Ting and Keller [6] 
(referred to as [TK]). The authors make a slender jet approximation where the radial length 
scale is taken to be much smaller than the axial length scale so that their ratio forms a small 
parameter which is used in an asymptotic theory to develop a system of one-dimensional 
coupled PDE’s. Similar ideas have also been used by Papageorgiou and Smith [7] in the 
study of the nonlinear evolution of plane wakes. The evolution equations of [TK] are simpler 
than the original system making them amenable to an analysis that produces families of 
pinching solutions which are physically relevant. Self-similar solutions in turn transform the 
PDE’s into ODE’s and numerical or exact solutions are given. 

The problem studied here is an extension of [TK]. In [TK] the cylindrical jet flows in air. 
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Within the framework of the Euler equations the undisturbed velocity field is a constant axial 
velocity which can be removed from the problem by a suitable Galilean transformation. The 
mechanism of three-dimensional Kelvin-Helmoltz instability is absent and surface tension 
dominates. Our interest lies in the competing effects of capillary instability and Kelvin- 
Helmholtz instability which is present if the surrounding (without loss of generality ambient) 
fluid is of non-zero density. Recently there has been a lot of interest in the formation of 
singularities in cylindrical vortex sheets. Linear stability studies include the work of Rotunno 
[8] and more recently singularity formation studies have been undertaken by Caflisch [9], 
Caflisch and Li [10], Caflisch, Li and Shelley [11], Grauer and Sideris [12], who do not 
include surface tension but include a background swirl in the azimuthal direction. The flows 
in [9]-[ll] do not possess a jump in axial velocity at the sheet surface and are shown in 
[8] and [11] to be linearly stable to short waves in the case of axisymmetric disturbances. 
The problem studied here omits swirl and includes axial velocity jumps which contribute a 
Kelvin-Helmholtz instability for axisymmetric disturbances. Axisymmetric flows are more 
likely to form pinching singularities since the axisymmetry prevents the core from flattening 
out. Even though we do not include swirl in our study we do enforce axisymmetry, and our 
slender jet equations predict self-similar pinching of the jet even in the absence of surface 
tension. Such solutions may be suitable candidates for initializing full-scale computations. 

The rest of the paper is organized as follows. Section 2 derives the slender jet equations. 
For completeness we derive a system for fluids of different densities with the axis of the jet 
being vertical. There, we also describe some useful transformations that reduce the num- 
ber of unknowns in the evolution equations from three to two which are also more relevant 
physically. Section 3 discusses properties of the equations including conserved quantities, 
equation type and linear stability characteristics. According to linear theory, stability win- 
dows exist for certain values of density ratios and certain undisturbed jet to pipe radius 
ratios. Equal density flows are always linearly unstable. In Section 4 we consider the possi- 
bility of finite-time singularities and derive the self-similar governing equations. In Section 5 
we solve these equations for vortex sheets in the absence of surface tension and give families 
of exact solutions, and in Section 6 we address the problem including surface tension. Section 
7 is devoted to discussion. 

2 Derivation of the evolution equations 

Consider the irrotational flow of an axisymmetric jet of density p\ and potential (f>i(t,i\z ) 
whose bounding surface is S(t,z). The flow outside is also irrotational and is taken to 
have potential <f> 2 (t,r,z) and density p 2 . The flow takes place in a vertical arrangement 
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with gravity acting and there is an outer bounding solid surface at r = /? 2 (a cylindrical 
polar coordinate system (r, 0, z ) is used throughout with the corresponding velocity field 
u = (u, u,u>)). If in addition surface tension acts and the tension coefficient is constant, a 
say, the following piecewise constant solution of the Euler equations is available: 

Ui = (0,0, Wo) 0 < r < Ri, 


u 2 = (0, 0, 0) R\ < r < Ri , 
cr 

Pi - Pi = 

where Wo, R\ are constants. This solution represents a perfectly cylindrical jet of constant 
axial velocity flowing int6 a quiescent medium. When the jet interface deforms we need to 
solve the problem in the regions T>\ : 0 < r < S{t,z) and P 2 : S(t,z) < r < Ri- There are 
two separate mechanisms of instability present: (i) Kelvin-Helmholtz instability of a cylin- 
drical vortex sheet, (ii) the destabilization of long waves by the circumferential component 
of the radius of curvature if surface tension acts. We describe the nonlinear evolution by 
assuming that the axial length scale of interfacial perturbations is large compared to the 
radial length scale. We note that such long length scales do not arise directly from linear 
stability analyses of jets (the characteristic length is of the order of the jet radius for capillary 
dominated flows), but enter the nonlinear evolution through initial interfacial conditions of 
the form f(tz) with e small. The following non-dimensionalizations are appropriate: 


r 



S l- Wo i 2 - 

% t ~~d t 


dW 0 


' 1,2 


(2.1) 


Utilization of tile scalings (2.1) into the governing equations (omitting the bars) gives: 
Laplace equations in V\^ 

<t>rr H <^r + ** <f> zz == 0 (^*2) 

T 

Kinematic conditions on r — S 


= £ 2 {St + <f>lzSz) 

hr = € 2 {S t + hzSz) 


Bernoulli equation on r — S 


<f>U + 2 (^2 ^lr + 4>\ z) ~ pi^it + y + $L) “ 2 ^° 


= -w e 


( 1 + e 3 ^)- 1 / 3 




s 


77 + 


i S 2 S~. 

rr M A,*, 


(1 + (1 + PSlfl* J ' 


(2.3a) 

(2.36) 


(2.4) 
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Boundary condition r = a 


<f> 2 r{r = a) = 0. 


(2.5) 


The constant terms in the Bernoulli equation (2.4) are chosen to satisfy the unperturbed flow. 
The term on the right hand side is the increment induced by surface tension on the pressure 
acting across the interface. In (2.2)-(2.5) above, the following nondimensional groups appear 


2 Pi ( 1 “Px 

E =1p e=rr(« = — ) 


p 1 


1 +P J 


W - - - — 

e R\ pi a R, 


These parameters are the ratio of radial to axial length scale, the density ratio between outer 
and inner fluid, the Weber number (ratio between capillary pressure and dynamic pressure), 
and the ratio between the cylinder radius and the unperturbed jet radius. Equations (2.2)- 
(2.5) are exact and are a result of a change of variables. For slender jets, however, e is a 
small parameter and the following asymptotic expansions are appropriate 


^ 1,2 = <^ 1°2 


+ e 


MS + 


(2.6a) 


S = So + ( 2 S\ + . . . . (2.66) 

Substitution of (2.6a,b) into (2.2) gives the following solutions that satisfy boundedness at 
the axis (?’ = 0) and the boundary condition at r — a 


4 0) = b,((,z) , 4 0) = &(m), 


(2.7) 


2 

= — ~^B\zz + C\(t, z) , 

= - j B izz + jB 2zz ln(r) + C 2 (t,z). (2.8) 

The functions B\ i2 and C\ t2 are to be determined from the kinematic and Bernoulli equations 
(2.3) and (2.4). The solutions (2.7) for the leading order potentials in T> \ <2 automatically 
satisfy (2.3a,b) and (2.4) to leading order. Next, at 0(c 2 ) and 0(1) respectively, equations 
(2.3a,b) and (2.4) along with the solutions (2.7), (2.8) give 


~^S 0 B Uz = Sot + B\ z Sozt 


(2.9) 


1 a 2 

— O SqB'Izz + ^ TrB2zz = Sot + B'lzSoz, 
l loo 


(2.10) 


1 \ w \ 

Bu + jB?, - KBit + 2 B'L) = + W.~ -Wl (2.1 1) 

Equations (2.9)-(2.11) as given above constitute a nonlinear system of three evolutions equa- 
tions for the three unknowns .S'o, B\, B 2 . Next we show how they can be reduced to a set of 
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two coupled equations for the shape of the jet and the axial velocity jump across it. This 
reduction is useful both mathematically and physically in discussing the pinching process 


after a finite time (see later). 

If we subtract (2.10) from (2.9) to eliminate S ot , tbe result can be readily integrated once 
with respect to z and enables expression of B i in terms of B 2 by the formula 


S*B U = (.S’o - a 2 )B 2z . 


( 2 . 12 ) 


The unknown function of t that arises from the integration can be taken to be identically 
equal to zero by choice of appropriate initial conditions (see Papageorgiou and Smith [7] also 
for a similar theory for two-dimensional wakes). Next we define new variables W and V 


given by 

W = B\ — B-i , V = B\ + B 2 , (2-1'^) 


which transform (2.12) into 

la'V, = (j - SI)W. . ( 2 - 14 ) 

The first evolution equation is found by addition of the kinematic conditions (2.9) and (2.10), 
and elimination of V while the second equation follows from (2.11) again after elimination 
of V using (2.14). It is convenient to differentiate the resulting Bernoulli equation with 
respect to z for two reasons: (i) The constant (which can also be a function of time in other 
configurations) associated with the Bernoulli integral disappears, (ii) the resulting equations 
involve the functions Sg and W z . The latter is seen to be the jump in axial velocity across 
the deforming vortex sheet and is therefore proportional to the sheet strength. Introducing 

the variables 2 

i , = % , X=W> , (2- 15 ) 

a 1 

gives the following coupled system of evolution equations: 

vt + {ix - n 2 x) z = ( 2 - 16 ) 


d + «(^ - v)]Xt + - r] 2 x). 

+ ( l f + l ( 5- ,)2 + 5 ( 5-’' )]x l = “^ 


(i + p) V 7 ? 


Equation (2.17) above, simplifies considerably for the case of equal densities a — 0; equation 
(2.16) is unaffected by density differences and (2.17) for a = 0 becomes 



Equations (2.16) and (2.17) (or (2.18)) are to be solved subject to prescribed conditions at 
t = 0. To fix matters, we will assume that the flow is periodic in z of period 2i r (a flow of 
any other period can be scaled to this canqnical case by rescaling time). 
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3 Properties of the evolution equations 

3.1 Conserved quantities 

The following are conserved quantities of the evolution 

r‘2n 1 l r'2* 

h= x(o + o a “ a7 /M~ > ^ = / rjdz. (3.1 a,b) 

Jo Z l Jo 

The first conserved quantity is a direct consequence of the Bernoulli equation (2.11). The 
total energy of the system is also a conserved quantity. The calculations that lead to (3.1c) 
below are included in Appendix A. 

h = J o (jp(l + -»l)V + ~V /3 ) dz. (3.1c) 

3.2 Equation type 

Equations (2.16) and (2.17) can be written in the form 


1 

X 


+ Mv 


4 l)r 0 ' 


where A is a 2 x 2 matrix with entries which are nonlinear functions of rj and x- It > s eas Y 
to show that the eigenvalues of A, A say, satisfy 


(X - %VX ~ X ) 


' (1 +oQx-2t/X 
2 A 


-*)- 


where 


T)(\ - Tj) (\ 2 2 We — 3/2 

—A— U (a - 1)x - anx -WT7) n 


* I ,1 
A = 2 +a( 2“’ ,) ' 


)■ 


The type of the evolution equations depends on the roots of the polynomial above. In what 
follows we show that the eigenvalues are complex, in almost all cases, and so the system is 
elliptic in general. To do this it is enough to show that the discriminant of the quadratic 
equation is negative. After a series of algebraic manipulations it is shown that the eigenvalues 
are 


where 


2A = x -2 TO + (1 + 0 ) X ~ 2TO *W 


D 2 _ 4xMl - q)( Q 2 - 1 ) 4W e 77(l - 77)77 3/2 


(1 + a — 2arj) 2 a(l -f p)(l + or — 2at]) 
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Since 0 < t? < 1 and a 2 < 1 we can conclude that the eigenvalues are complex since D 2 < 0. 
The reason for this is that both terms in the expression for D 2 above are negative; this 
conclusion uses the fact that 2 A = (1 + a — 2 arj) > 0 for all relevant values of a and rj. 

It is concluded, therefore, that the system is elliptic for almost all values of cv and W e . 
The exceptional case has a = 1, W e = 0 which represents a jet in air with zero surface tension 
yielding real but repeated eigenvalues (further discussion of this case is given in Section 7). 

The linear stability spectrum follows from a linearization of the system (3.2) about the 
unperturbed states. If we write 

(x,i) = (U^) + s(x,vV H ’- a \ 

it follows that c is given by ±(Z) 2 ) 1 / 2 where D is as given by (3.3) and evaluated at (x, 7) = 
(l,4f). As shown above, D 2 < 0 and so instability is obtained with growth rates of 0(A;), 
as in Kelvin-Helmholtz instability of a flat vortex sheet for instance. Our aim here is to 
follow this instability into a fully nonlinear regime with the cylindrical jet interface growing 
to amplitudes comparable to the unperturbed radius. 


4 Finite-time singularities: Jet pinching 


In this section we consider analytically possible nonlinear terminal states of the evolution 
equations. It is found that density differences are unimportant in the leading order balances 
between unsteady and nonlinear terms as a singularity is encountered. It is sufficient to 
consider equal densities p = 1 (a = 0). The case p 7^ 1 also leads to equations (4.4a,b) 
below, with W e replaced by 2 W e . For equations (2.16) and (2.18), therefore, we postulate 
that the solutions encounter singularities after a finite time, t s say, at axial positions z s and 
seek to describe the solutions near such points. We make the following ansatz 


rj(z, t ) = (t, - *)“/(0 x(*> t) = (t s - <)M0 

( z ~ z ») (4.1) 

? (ts-tf ^ ’ 

Physically relevant solutions have « and ft positive in order to ensure that the interface 
position does not become unbounded at the critical time, and that the dominant features 
are due to nonlinear focussing. Under the transformations (4.1) time and space derivatives 
become 

fit 8 d_ ^ 1 d 

(ts ~ t ) d£ dz 


d 8 

y L 

dt dt 


(4.2) 


(ts-tydf 

Substitution of (4.1), (4.2) in (2.16) and (2.18) along with the usual order-of-magnitude 
arguments as t —* t 3 from below gives the following expressions relating the constants: 

a = 4 — 4/3 7 = (3 — 1, (4-3) 
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and for a to remain positive we require 


0 < (3 < 1. 

Note that this choice of (3 implies a negative 7 and so the axial velocity jump blows up as 
the jet radius shrinks to zero. 

The conserved integrals (3.1a-c) become to leading order in ( t„ — t ): 

/ oo - /*oo 

9 d( h = ( 1 , - ty-w fd( 

-00 J -00 

h = (i. - t)‘- p r ^ ( 4 ' 4 > 

J —00 4 (l 

Equations (4.4) above represent the local contributions of the focussing region to the three 
conserved quantities. If the exponents of ( t s - t) are positive the contributions are zero 
in the limit and no integral constraints need to be placed on possible solutions f(() and 
( /(£), whereas if any exponent is negative or zero the corresponding integral it multiplies 
must vanish. Using (4.3a) we see that the I 2 and I 3 contributions are zero in the limit. 
Consideration of the h contribution, however, shows that the following constraint is required: 

/ OO 

<7^ = 0. (4-5) 

-OO 

The governing equations transform to the following one-parameter system of nonlinear 
ODE’s: 

-4(l-W + (W + (f9)' = 0 ’ ( 4 - 6a ) 

(1 - (3)g + f3(g' + l -{g 2 )' = ~(r' /2 )' . (4.66) 

In general (4.6a,b) must be solved numerically for different values of /? and with conditions 
at £ = 0 (i.e. z = z s ) specified for example. 

Before the presentation of solutions to (4.6a,b), we note that the same self-similar equa- 
tions are obtained for the case of a jet in air. Further more, the constraints on (3 imposed 
by the conserved integrals of the system are preserved and so the solutions that follow ap- 
ply equally well to jets in air. For completeness the analogous derivation of a jet in air is 
contained in Appendix B. 

5 Collapse of cylindrical vortex sheets with zero sur- 
face tension. 

The equations governing the self-similar break-up structure become: 

-4(1 - /J)/ + «/' + (/«)' = 0, (5.1a) 
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(1 - 0)9 + M + 5(s 2 )' = 0- (' r >'6) 

The equation for g is nonlinear but is decoupled from that for / which is linear. Once g is 
determined, the solution for / readily follows: 


fit) = fo 



5$ + 4# 


m + a) 




(5.2) 


The main task, therefore, is in the determination of g and classes of exact solutions are 
constructed next. 

An exact solution is readily available for the critical case ft = \', (5.1b) can be integrated 
once to yield a quadratic algebraic equation for g(£) which yields the solutions 


2<7± = — £ ± \A 2 + 4<7o> 


(5.3) 


where g 0 = g(0). The solution <7+ pertains to positive velocity jumps while g. describes 
flows with negative jumps. The corresponding solutions for / are 

fo (( + yfr + ^ 

U 16j;o « 2 + *9l)' n 

f- = 64/o0o 6 (« + l/C 2 + 4s 2 o)- 5 (? + 4 (5.4) 
with fo = /( 0). The behavior for large £ is 

r fo t4 9o £ v 1 „ 

/+ ~ T7 - 4^ <7+ ~ T * ^ +°°> 

16^o £ 

/+~2 5 ®|^- 6 g+ ~ |£| £->-oo. 

/- ~ 64^q/o^ -6 <7- ' 1 £ -» +°°, 

5o ki 

Even though the solutions /±(£) are admissible the corresponding <7±(<0 are not, due to 
the conservation constraint (4.5) requiring that the integral vanish. This is satisfied by 
</(0 =5-(0 for £ < 0, and g(£) = g+(() for £ > 0, but such solutions are not acceptable since 
g(£) and /(£) are now discontinuous at £ = 0, in violation of our slender jet approximation. 
We reiterate here, that any solutions /(£), g(£) with bounded derivatives are acceptable 
singular states of the slender jet equations; unbounded derivatives herald the introduction 
of shorter scales over which a new acceptable solution may emerge. Such analyses are not 
pursued further here since we are interested in the construction of pinching solutions. 

For any value of ft, 

««) = -< no = 
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are exact solutions with A some constant. The solution for g is an odd function of £ and is 
consistent with the integral constraint (4.5). For /(£) to be acceptable we require bounded- 
ness of the first derivative at £ = 0; this is not possible for ft < 1, and this family of solutions 
is therefore dropped. 

For ft — 2/3 another exact solution is found; this is most easily seen by multiplication of 
(5.1b) by g and integration with respect to £ which gives 

(i - |/J) ft + \k<j‘ + r/ = jsl so = s(0), 

from which it follows that for ft = 2/3, g satisfies the cubic polynomial 

i + 4s 2 - d = 0. (5.5) 

If g 0 — o the nontrivial solution is g = —( with corresponding / = A which is singular 
at the origin and is discarded (see above). Equation (5.5) can have (i) one real root and a 
pair of complex conjugate roots, (ii) three real roots of which at least two are equal, and (iii) 
three distinct real roots, depending on the conditions 

(i) Q > 0 (ii) Q = 0 (iii) Q < 0, 

Q = b/o(<7o - ), 

respectively. The position £ = £<- = ^ 3 < 7 o yields Q — 0 and is the point above which three 
distinct solutions exist. At ^ there is a root g = <7o/4 lyr3 and a pair of equal roots 

g = — 2</o/ 4 1/,:3 . A local analysis shows that g ~ (£ — l) 1 ^ 2 in the vicinity of the double root 
indicating that the derivative is infinite at . The solutions of (5.5) for g$ — 4/27 are shown 
in Figure 1; the dashed line represents g = — £ which is the asymptotic state of one of the 
solutions for large £. Within the context of the long wave theory, the roots with unbounded 
derivatives are discarded and so a unique acceptable solution exists for ft = 2/3. 

The type of solution given above for ft = 2/3 can be generalized to arbitrary values of ft. 


Under the transformation 

g (0 = (G(o, 

(5.6) 

equation (5.1b) becomes 

T r 



Gp(\ + Gy - 0 = j, 

(5.7) 


where K is a constant of integration. Equation (5.7) yields G implicitly in terms of /. We 
can rewrite (5.7) in terms of g using (5.6) to give 



where go = g{0). Equation (5.8) is a generalization of (5.5) to arbitrary values of ft. In 
particular if ft is a rational number with ft = p/q with p, q integers such that 1/2 < ft < 1, 
(5.8) becomes 

S p (( + g)'-” = 9l ( 5 - 9 ) 

The values p = 2, (7 = 3 yield equation (5.5). We are interested in generating real solutions 
of (5.9) which are acceptable similarity forms of the breakdown. To do this we begin at 
£ = 0 where there is one real root and q — 1 complex roots if q is odd, or, two equal and 
opposite real roots and q — 2 complex roots if q is even, and follow the appropriate real root 
as a function of £ to obtain <jr(£). This can be done for negative £ also; alternatively the 
differential equations can be solved directly. This was also done as a check. 

The behavior of solutions to (5.1a,b) for large £ is found to be 

/ ~ (-+oo. (5.10a) 

The constants of proportionality can be determined numerically. For £ — ► — oo the equations 
are analyzed for large positive £1 = — £ to give the asymptotic forms 

/~(- 0 " 1 9 =*{-() + •- Z^-oo. (5.106) 

We note that the constant of proportionality for the decay rate of / is determinable numeri- 
cally, whereas the constant for the leading order asymptotic behavior for g is equal to unity. 
These results are confirmed by our numerical solutions described next. 

Numerical solutions can be generated for any ft in the designated range. Here we present 
a set of solutions for ft = 2/3,4/5,9/10. The results are presented in Figures 2, 3 for /(£) 
and g(£) respectively. The initial conditions for all solutions are /( 0) = 0.5, g(0) = 1.0. 
The asymptotic results (5.10a,b) are confirmed by these solutions. The behavior of g(£) as 
£ — ► — oo predicted above is clearly seen in Figure 8 where the behavior is independent of 
ft and the curves collapse to a line of unit slope. Of particular interest are the solutions for 
ft = 4/5. According to (5.10a), /(£) grows linearly for large positive £ when ft — 4/5. For 
ft < 4/5 the /(£) is concave for large £ and for ft > 4/5 it is convex. As explained below, 
viewed from outside the similarity region, the solution for ft = 4/5 represents a conical jet 
(the radius is a linear function of axial position) while the solutions for ft < 4/5 and ft > 4/5 
represent concave and convex jets respectively at the singular time. The convex jet solution 
which has ft = 7/8 is shown in Figure 4 for £ > 0. The initial conditions are /( 0) = 0.3, 
<j(0) = 0.6. The asymptotic forms (5.10a) are confirmed numerically in Figure 5 by use of 
a logarithmic plot and a least squares fit of the resulting data to obtain an estimate for 
the slope. As expected the errors between the exact and the computed value of the slope 
decrease slowly as the upper limit of integration is increased. 
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Since the solution for / becomes unbounded for large positive £ (as does that for g at 
large negative £), it is essential to check that these asymptotic values are consistent with 
the ansatz (4.1); such analysis also provides a qualitative description of the break-up from 
outside the similarity region. If we substitute the asymptotic forms (5.10a) into (4.1) and 
express £ in terms of the physical spatial variable z — z s (> 0) we find that when viewed from 
outside the self-similar region, the break-up takes the spatial characteristics 

4(1— 0 ) 1 — J8 

T) ~ (Z — Z a ) P X ~ { z ~ z s)~ 0 ■ 

For z < z s the asymptotic forms (5.10b) are used to give 

- t) 4 -^ (z t - z) 

Z 3 - Z) X (t s ~ t ) 

These expressions indicate the geometry of the jet at break-up; the shape is squeezed to a 
thread as £ — > — oo with the radius tending to zero as the singular time is approached (this 
can be seen from the fact that 1 < 4 — 3/? < 5/2) on one side, while it has a concave or 
convex growth on the right depending on the value of ft. The corresponding velocity jumps 
tend to infinity and zero respectively at the algebraic rates given above. For example, the 
exact solution given above for ft = 2/3 gives a parabolic shape at break-up. 

6 Collapse of cylindrical vortex sheets including sur- 
face tension 



The governing equations are (4.6a,b). The parameter W e /a in (4.6b) can be re-scaled to 
unity by the transformation / — * (kk 2 /a 2 )/, yielding the normalized form 


-4(l-0)/ + /^/' + (M = O, 


(6.1a) 

( 6 . 16 ) 


Before presentation of numerical solutions it is useful to examine the large [£| behavior of 
the solutions. For large positive £ we look for solutions of the form 


/ ~ v4C l 

Using these expressions in (6.1a) implies 


g ~ B( " m, n > 0. 


( 6 . 2 ) 


m = 


4 ( 1 -/?) 
ft ' 
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The correct value of n which matches the numerical solutions is found by a balance of the 
linear terms in (6.1b). To illustrate this we write the asymptotic form of (6.1b) according to 
the scalings (6.2), 

(i _ - pnBc n - B 2 nc 2n ~' + o(r 2n_1 ) = 

with part of the higher order correction terms shown explicitly. If we choose to match the 
linear terms in B on the left with the surface tension nonlinearity on the right, a value for 
n is found which makes the constant B negative in contrast to numerical solutions. The 
correct behavior is determined by matching the linear terms which gives 


m 1 — (5 



With this value of n the linear terms (which are the largest terms in the equation) match 
automatically with B arbitrary. Similar arguments apply for £ — > — oo (the analysis is 
almost identical if we write and analyze the equations as — ► oo). To summarize, 

therefore, we find 

g ~ Z -* °° (6.3a) 

/-(-O 1 ^ t^-oo (6.36) 

Note that the solutions for large positive £ have the same asymptotic behavior as those corre- 
sponding to zero surface tension. As £ — » — oo, however, the behavior is altered significantly. 
The consequences of this are discussed after the presentation of numerical solutions. 

Equations (6.1a,b) were solved numerically with particular interest in the verification of 
the asymptotic structures (6.3a,b). The initial conditions used were /( 0) = 1/2, g(0) = 
1 with an integration range as large as ( — 100, 100). Results are presented for the three 
representative values ft = 2/3,4/5,9/10. Figure 6 shows collectively the solutions for g(£). 
The maximum/minimum in g decreases/increases as ft is increased, A notable feature is 
that there are two distinct regions, £ < £o, £ > Co f or some value £o, where the axial velocity 
jump is negative and positive respectively. Changing the initial condition so that g{0) < 0 
for example, indicates that this separation is preserved. Results for the corresponding / 
solutions are given in Figures 7a,b,c separately due to the large scale differences involved. 
For ft < 4/5 the computed shape is concave, and the solution for ft = 2/3 has a minimum at 
£ = —0.3 approximately; the solution corresponding to ft = 4/5 is a critical one with a linear 
growth at ±oc (see Section 5 also), with a minimum radius at £ = -0.38 approximately; 
solutions having f3 > 4/5 are convex in regions sufficiently far away from the origin and 
the computed solution for ft = 9/10 has a minimum at £ = -0.48 approximately. For the 
solutions presented here, the position of minimum radius (i.e. the value of £ which makes 
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/(/) a global minimum) is larger than the positions where the jump in axial velocity (which 
proportional to </(£)) is zero - 

Next we use the numerical solutions just presented to check the asymptotic forms of the 
solutions. This is done using logarithmic plots and a least squares fitting of the resulting 
data. For large positive £ results are shown in Figures 8a-c for / and g for the three 
representative values of ft. The agreement is very good and is found to improve as the value 
of computational “infinity” is increased, as expected. For large negative £ analogous results 
with similar accuracy are found and need not be presented here. 

As in Section 5, the asymptotic form of the solutions far from the origin enables the 
qualitative description of the shape at break-down. Considering the shape first, there appear 
to be three canonical forms at break-down corresponding to the cases (i) ft < 4/5, (li) 
ft = 4/5, and (iii) ft > 4/5. For case (i) a concave shape is found at infinity and the jet 
would tend to pinch forming drops; case (ii) gives linear shapes at infinity and at break- 
down the jet surface would look like two right circular cones with a common vertex; case 
(iii) yields jets with shapes which are locally convex in the axial direction. Photographs of 
experiments in air (see [4], [5]) indicate behavior which fits in with cases (i) and (iii). An 
additional fact which is evident from the experiments and makes long wave theories like the 
present one (see also [TK]) useful candidates for break-down prediction, is the validity of the 
slender jet assumption and in particular the absence of gradient steepening in the shapes at 
break-down. 

7 Conclusions 

We have developed a nonlinear theory to describe the pinching process of a cylindrical jet 
of fluid surrounded by a second fluid with possibly different density. The interface between 
the two fluids can support surface tension. The theory is based on the smallness of the 
jet radius to a typical axial wavelength of the jet. This excludes jets which try to roll up 
and such behavior would appear as a slope singularity in our model evolution equations. 
As an illustration of such behavior, consider the equations (Bl), (B‘2) for a jet in air in 
the absence of surface tension ( W e = 0). Equation (B2) is an inviscid Burgers equation 
and x forms a shock singularity after a finite time (for definiteness we are considering the 
periodic problem with smooth initial conditions). Similar behavior is found for 17 ; the slender 
jet approximation is violated for such solutions since the axial length scale becomes much 
smaller than the radial one violating the original assumptions used to derive the model. If 
surface tension is included, however, the equations become elliptic and allow the possibility 
of admissible pinching solutions described in this article. It is interesting to find that for jets 
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surrounded by an outer fluid so that the Kelvin-Helmholtz instability mechanism is present, 
the model equations even in the absence of surface tension are elliptic and pinching solutions 
can again be constructed. 

Finally, we return to the similarity solutions constructed in Sections 5 and 6. There 
appear to be three canonical breakdown regimes corresponding to (i) ft < 4/5, (ii) ft = 4/5, 
and ft > 4/5 (see Section 6). In what follows we consider the validity of such solutions within 
the context of the long wave theory. These equations are derived for the following unsealed 
lengths (stars are used to denote unsealed quantities): 

z* ~ — r* ~ 1 . (7. 1 ) 

e 

As a finite-time singularity is encountered, it is seen from (4.1) that the length scales over 
which similarity solutions are constructed are: 

„ ft* ~ ft . r * ~ (t s - t) a . (7.2) 

We know from (4.3) that a = 4 — Aft which when substituted into (7.2) along with the 
ordering z* » r* in order to validate the slender jet approximation, yields the ordering 

» ( t s -t) 4 - 4( * . (7.3) 

The three canonical cases (i)-(iii) described above can be seen in terms of (7.3). It is eas- 
ily verified that cases (i) and (ii) satisfy (7.3) for all t arbitrarily close to t s . Case (iii), 
characterized by ft > 4/5, is analyzed by letting ft = 4/5 + 8 with 8 > 0; this requires the 
ordering 

(t s — t) 5S >> e or (t s — t) » c*?. (7-4) 

Equation (7.4), therefore, defines a range of times for which convex interfacial shapes are 
formed at break-down. From a physical point of view we expect solutions corresponding to 
cases (i) and (ii) to be more appropriate when the ordering (7.4) is not observed. Composite 
or patched solutions are possible, however, whereby for times given by (7.4) a convex drop 
begins to form which is then matched to a concave or linear shape corresponding to ft < 
4/5 and ft = 4/5 respectively, to allow evolution well into the break-down regime without 
violation of the slender jet ansatz. 
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A Derivation of the energy conserved quantity 

The flow is governed by the incompressible Euler equations which are conservative and 
so the total energy of the system is expected to be conserved throughout the evolution. 
Contributions to the energy come from (i) kinetic energy, (ii) potential energy (in the presence 
of gravity), (iii) interfacial energy. For a periodic flow of period d (non-periodic, flows simply 
alter the integration limits) the total energy using the dimensional variables of Section 2 is 


r R 2 fd/2 1 fd/2 1 - .0 , 7 

S= I II -p 2 \V<h\ rdrd0dz+ / / / -pi\V(f>i\ rdrdOdz 
jS Jo J —d/2 2 JO Jo J—d/2 £ 

+ J a I L^ rzdrd)dz+ i L L 


fd/2 

-\-(T / / Sdddz. 

Jo J —d /2 


(Al) 


Clearly the potential energy contributions are zero by virtue of the choice of reference frame. 
Non-dimensionalizing (Al) according to (2.1) (and dropping bars) yields 

*-fCMM'***ra»*«)’** 


r 1/2 

+w e / ; 

J- 1/2 


S dz. 


(A 2 ) 


Making use of expansions (2.6a,b) along with the leading order solutions (2.7) and the change 
of variables (2.15) yields the following conserved quantity 


So = f'^ (lp(l - v)v 2 x 2 + v) 2 x 2 + ^V /2 ) di 


(A3) 


where £ 0 is a constant. The limits of integration are normalized to be 0 and 2 t r as they 

appear in equation (3.2). ... .. . _ 

We have shown, but for brevity do not include here, that (A3) is consistent with the 
evolution equations (2.16), (2.17) for the equal density case. This is achieved by taking 
the time derivative of (A3) and eliminating t] t and Xt from (2.16), (2.17) to show that the 
integrand is an exact z-derivative and so the expression is identically zero. 
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B Equations for the collapse of a jet in air 

The analysis of Section 2 is used with a zero density surrounding fluid which implies that 
(j ) 2 = 0. The equations to be solved, to leading order, are (2.9) and (2.11) with - 0. If 
we define 

V = S 2 0 X = 


the coupled system becomes 

Vt + {vx)z = 0, {B\) 

X, + = -W,(r' r % (B2) 

The equations are in conservative form and so the integrals (over space) of rj and \ are 
conserved quantities (c.f. (3.1a,b)). The energy conserved quantity is easily obtained from 
the results of Appendix A and in particular equation (A2). The result is 

£, = Jo* (i 7 ** 2 + We7?1/2 ) dz - (53) 

We note that (Bl) and (B2) are equivalent to equations (2.11) and (2.10) of [TK] respectively 
if (2.10) is differentiated with respect to z and the equations are put in terms of our variables. 

As a check, if the time derivative of (B3) is taken and rj t , Xt are eliminated by use of 
(Bl), (B2), the resulting integral is 

-t / [(v X 3 )z + W e (ri 1/2 x)z]dz, 

4 Jo 

which is zero, as expected, since it is an exact derivative. 

If we look for singular solutions of the form (4.1), the orderings (4.3) emerge as well 
as the scaled conserved quantities (4.4). Mathematically, therefore, the two problems are 
equivalent and the same restrictions on (J apply. The equations governing collapse in the 
self-similar region are, to leading order, the system (4.6a,b). 
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Vortex sheets zero surface tension 



FIGURE 1 Similarity solutions. Zero surface tension, W e = 0. Possible solutions for 
scaled axial velocity jump g(£) for (3 = 2/3. 
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FIGURE 2 Similarity solutions. Zero surface tension, W e = 0. Scaled interfacial shape 


/(£) for different (3 — 2/3,4/5,9/10. 
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We=0 log(f) , log(g) 



log(g) - siope»0.56934 (exact 0.571429) 


log(f) - slope--0.142142 (exact -0.142857) 


‘ 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9 

log(xi) 

FIGURE 5 Similarity solutions. Zero surface tension, W e = 0. Logarithmic plots to 
show the asymptotic behavior of /(£) and g(£) for large £ and fi = 7/8. Slopes computed 
by least squares included on the figure. 

We/a=1 g for different beta 


FIGURE 6 Similarity solutions. Non-zero surface tension, W e - 1. Scaled axial velocity 
jump gtt) for different /? = 2/3, 4/5, 9/1 (K 







FIGURE 7 Similarity solutions. Non-zero surface tension, W e = 1. Scaled interfacial 
shape /(£) for different ; (a) /3 = 2/3, (b) /3 = 4/5, (c) /? — 9/10. 









We/a=*1 log plots beta=2/3 

xl -> + infinity 
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We/a=1 log plots beta=4/5 


xi -> + infinity 




FIGURE 8 Similarity solutions. Non-zero surface tension, W e = 1. Logarithmic plots to 
confirm asymptotic, states for large £. Behavior of /(£) and g(£) for different ft; (a) ft = 2/3, 
(b) ft = 4/5, (c) ft = 9/10. Slopes computed by least squares included on the figure. 
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